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' Context. The onset of gravitational collapse in cluster-forming clouds is still poorly known. 

, Aims. Our goal is to use the Class protostar IRAS 4A, which is undergoing collapse in the active molecular cloud NGC 1333, to set constraints 
on this process. In particular we want to measure the mass infall rate and investigate whether the collapse could have been triggered by a strong 
external perturbation. 

Methods. We analyze existing continuum observations to derive the density structure of the envelope, and use our new molecular line observa- 
Q ' tions done with the IRAM 30m telescope to probe its velocity structure. We perform a detailed comparison of this set of data with a numerical 
^ model of collapse triggered by a fast external compression. 

CZ3 Results. Both the density and velocity structures of the envelope can be well fitted by this model of collapse induced by a fast external compres- 
' sion for a time elapsed since point mass formation of 1-2 x 10'' yr. We deduce a large mass infall rate of 0.7-2 xlO""* Mq yr"' . The momentum 
^ , required for the perturbation to produce this large mass infall rate is of the same order as the momenta measured for the NGC 1333 numerous 
? ""j ' outflows. Our analysis shows also that the turbulence is highly non uniform in the envelope, dropping from supersonic to subsonic values 
l^sj ' toward the center. The inner subsonic turbulence is most likely a relic of the conditions prevailing in the dense core before the onset of collapse. 
Conclusions. The vigorous collapse undergone by IRAS 4A was triggered by a fast external compression, probably related to the expansion of 
a nearby cavity, which could have triggered the collapse of the nearby Class protostar IRAS 4B simultaneously. This cavity could have been 
generated by an outflow but we have not found a good protostellar candidate yet. 
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1. Introduction 

While the first phases of protostellar collapse in distributed 
star forming regions like the Taurus molecular cloud start to 
be observationally better constraine d (e.g. Tafalla et al. 1998; 
lOnishi et alJl999(lMotte & Andr e 2001; B elloche et alJ2 002). 
little is still known about the onset of gravitational collapse in 
cluster-for ming clouds where this process is likely to be more 
violent (see lAndre et alJ2004l) . 

The protostar IRAS 4A is located in the active molec- 
ular cloud NGC 1333 forming low- and intermediate-mass 
stars in the Perseus complex (e.g. Isandell & Kne3 [200 ll) . 
The NGC 1333 region contains a double cluster of infrared 
sources which is one of the least evolved embedded clus- 
ter known so far, with an age of 1-2 xlO* yr jLada et al] 
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1199(1 lLada&Ladall2003h . It contains also a population of 
younger stellar objects, including a few Class protostars (e.g. 
Looney. Mundv. & Welc Hl2000l) . IRAS 4A, one of these Class 
protostars ( Andre et alJll993h . ha rbors a 1.8" Ipinary sys- 
tem at a position angle RA. of -50° jReipurth et alJl2002l) and 
is located in the vicinity of another yo ung multiple system, 
IRAS 4B ( Lav. Carlstrom. & Hilljll995 '). Both sources are as- 
sociated with molecular outflows oriented approximately i n 
the North-South direction at 10"-scale (e.g. lChoilbOOlLbOOSt) . 
On larger scale, t he IRAS 4A outflow has a RA. - 45° (e.g. 
lBlakeetalJfl995l) . The whole region is actually filled with 
about ten molecular outflows driven by you ng protostars be- 
longing to the NGC 1333 protocluster (e.g. iKnee & Sandelll 

2m)- 

IRAS 4A was i dentified as a goo d infall candidate in th e 
surveys of .Mardones et al.l (1 19971) and lGresersen et all (Il997l) . 
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iDi Francesco et alJ ( l200ll) used the IRAM Plateau de Bure in- 
terferometer to probe the inner parts of its envelope. They de- 
tected inverse P-Cygni profiles in H2CO(3i2-2n) and CS(3-2) 
which they interpreted as infall motions. They derived a large 
mass infall rate of 1.1 x 10""* yr"', which is about 70 
times larger than the standard accretion rate ^ at 10 K jShd 
Il977l) . Such a large mass infall rate, if confirmed, can not oc- 
cur in an envelope collapsing spontaneously. The collapse of 
the IRAS 4A envelope was therefore very likely triggered by 
a strong external perturbation. Several authors already argued 
that the numerous outflows have created cavities affecting the 
density structure of the molecular cloud and eventuall y trig- 
gered further star forrnation ( Warin et al 1996: Lefloch et alJ 

ll998tlKnee&Sandel|E00(i lOuillenetal.ll2005l) . In this re- 
spect. ISandeU & K nee' ("200 1*) proposed that NGC 1333 is "an 
example of self-regulated star formation". 

However, it still remains to be shown that a collapse trig- 
gered by an external perturbation can match the density and 
velocity structure of a protostellar envelope such as IRAS 4A. 
This is the purpose of this work. The layout of the paper is as 
follows. Sect. 121 summarizes observational details. In Sect. |3] 
we derive the density structure of the envelope from a com- 
pilation of existing continuum data and we interpret in terms 
of velocity structure our new molecular line observations done 
with the IRAM 30m telescope. We compare these results in 
Sect. 0]to radiative transfer models of hydrodynamical simu- 
lations of collapse triggered by a fast increase of the external 
pressure ( He nnebelle e t al. 2003, 2004). Finally we discuss the 
implications in terms of triggered star formation in Sect. |5] 

2. Observations 

We carried out milimeter line observations with the IRAM 
30m telescope at Pico Veleta, Spain, in September and October 
2001, and in August 2004, in the following molecular tran- 
sitions: HCO+(1-0), Hi3CO+(l-0), HC"*O+(l-0), N2H+(l-0), 
CS(2-1), 0348(2-1), HCN(l-O) at 3 mm, CS(3-2), C3'*S(3-2), 
N2D+(2-l) at 2 mm, H2CO(3i2-2n), Hi3CO(3i2-2„), CS(5- 
4), N2D+(3-2) at 1.3 mm, and HCO+(3-2), Hi3CO+(3-2) at 

1 mm. The references to the frequencies we used are given in 
Sect. 13.21 The half -power beamwidths can be computed with 
the equation HPBW (") = We used four SIS heterodyne 
receivers simultaneously and an autocorrelation spectrometer 
as backend. The spectral resolution was 20 kHz at 3 mm and 

2 mm, and 40 kHz at 1.3 mm and 1 mm in 2001, and 10 kHz at 

3 mm, and 20 kHz at 1.3 mm in 2004. The observations were 
done in single-sideband mode with sideband rejections of 0.01 
at 3 mm and 0.05 at 2, 1.3 and 1 mm. Accordingly, the cali- 
bration uncertainty was ~ 10%. The forward efficiencies F^^ff 
were 0.95 at 3 mm, 0.93 at 2 mm, 0.91 at 1.3 mm and 0.88 
at 1 mm. The main beam efficiencies were computed using the 
Ruze function B^b = 1.2e e-(4'f«'^/'i)', with e = 0.69, Ro- = 0.07 
and A the wavelength in mm. The system temperatures ranged 
from ~ 110 K to ~ 160 K at 3 mm, ~ 190 K to ~ 250 K 
at 2 mm, ~ 320 K to ~ 670 K at 1.3 mm and ~ 550 K to 
~ 1200 K at 1 mm. The telescope pointing was checked every 
~ 2 hours on Saturn, 3C84, and/or NRAO 140, and found to 
be accurate to 3 " (rms). The telescope focus was optimized 



Table 1. Interferometric measurements of dust continuum 
emission. 



Ref." 


» A 


HPBW 


peak 




Size 


Interf. 




(mm) 


(" X ") 


(mJy.beam"^ ) 


(mJy) 


(" X ") 




DOl 


1.33 


2.0 X 1.7 


1200 ± 16 


3100 ±470 


5.3 X 4.1 


PdBI 


DOl 


2.04 


3.2 X 2.3 


450 ± 14 


690 ± 110 


4.4 X 3.4 


NMA 


DOl 


2.19 


3.4 X 2.6 


420 + 10 


580 ± 90 


4.2 X 3.3 


NMA 


LOO 


2.7 


5.5 X 5.0 


351 ±3 


544 ± 14 


25 X 24 


BIMA 






3.0 X 2.8 


280 ±2 


526 ±9 


12 X 18.5 


BIMA 






1.2 X 1.1 


172 ±2 


450 ± 10 


5.4 X 6.2 


BIMA 


DOl 


3.22 


2.9 X 2.6 


200 ± 4.4 


320 ± 49 


4.5 X 3.7 


PdBI 


DOl 


6.92 


2.1 X 1.9 


16 ± 0.29 


26 ±4 


3.3 X 2.8 


VLA 


DOl 


13.3 


3.8x3.1 


3.1 ±0.27 


4.7 ± 1.0 


6.3 X 3.6 


VLA 



Refere nces: LOO: iLoonev et alj llOodt) : DOl: iDi Francesco et all 
j200lh . 

Quoted errors are statistical errors only and do not include any 

systematic calibration uncertainties. 

Flux integrated over a region of size listed in column 6. 



on Saturn and/or 3C84 every ~ 3 hours. Position-switching 
and on-the-fly observations were done with a reference po- 
sition located at either (3000",3000"), (1500",1500"), (0",- 
240"), or (150", 120") relative to the envelope center q;j2ooo - 
03h29m^0.^30, 6o( m ^ 31° 13^3 1 .8", as measured in the 1.3 
mm emission by Motte & Andrei i2001j '- The data were re - 
duced with the CLASS software package ( iBuisson et alJ20oA . 
The spectra were converted from antenna temperature to main 



t:. 



the JCMT archive 
Hi3CO+(3-2) and 



beam temperature using the equation T^b 
In addition, we downloaded from 
CS(5-4), CS(7-6), C34s(5.4), HCO+(4-3) 
DCO^(5-4) spectra observed toward IRAS 4A in February and 
July 1992, February 1994, November 1999, and Octo ber 2001 
The data were conve rted into fits files with SPECX (Pad manI 
.1993: Matthews (feJenness, . 1997.) . reduced with CLASS and 
converted to main beam temperatures using beam efficiencies 
7;^b = 0.69 at 1 mm and 0.63 at 0.85 mm (see lMaret etalJ2004l 
for the conversion to Jmb at JCMT). 

3. Results 

3.1. Continuum emission 

The mass distribution of the IRAS 4A envelope can be es- 
timated from the dust continuum emission. We compiled 
from the literature the interferometric (Tab. Q and single-dish 
(Tab. 12} measurements done so far Masses were estimated 
by various authors using different assumptions for the dust 
emissivity, the dust temperature distribution and/or the dis- 
tance to NGC 1333. Here, we analyze this set of data ho- 
mogeneously. We assume a distance of 318 pc, based on the 
Hipparcos parallactic measurements of the Perseus OB2 a sso- 
ciation (Ee Zeeuw et al. 1999), as discussed by Getman e t alJ 
ll2002h . With this distance, the bolometric luminosity derived 



' Note th at our central position is offset from the two components 
observed bv lReipurth et alJ l2002h with VLA by 3" and 1.5" toward 
the West. 
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Fig. 1. Mass distribution estimated from the integrated flux 
of continuum emission measured with a) interferometers (see 
Tab.^ and b) single-dish telescopes (see Tab.|2j, and computed 
using the best-fit opacity exponent yS = 0.85. 

Table 2. Single-dish measurements of dust continuum emis- 
sion. 



(1) 



(2) 



(3) 



Refe rences: L98: lLefloc h et alJ tl998): SOO: ISmith et alJ f20QO): 
S01: ISandell & K nee 1 20dll): MOli lMotte & Andrd <200lh . 
Quoted errors are statistical errors only and do not include any 
systematic calibration uncertainties. 
Flux integrated over a region of size listed in column 6. 



this approximation should be valid, provided the density pro- 
file follows a power law. For larger radii, we expect the dust 
grains to be also heated by the external UV radiation field via 
the photoelectric eff'ect (e.g. Lesaffre et al. 2005). For simplic- 
ity, we assume that the dust temperature stays at 10 K for a 
radius larger than 9200 AU. In this case, the correction factor 
to compute the mass within a radius R > 9200 AU is 



■p-q 



(I) 



3-p \-l 
+ 1 



(2) 



Ref.<" A 

(mm) 


HPBW 

(") 


(Jy.beam"') 


int 

(Jy) 


Size 

(" X ") 


Telescope 


SOO 0.85 


16.0 


10.30 ± 0.03 


12.9 


30 X 23 


SCUBA 


SOI 0.85 


14.0 




9.05 


14.2 X 5.4 


SCUBA 


L98 1.25 


11.0 


3.10 ±0.01 


4.7 


18.8 X 18.8 


30m 


MOl 1.25 


11.5 


4.10 ±0.02 


4.1 


24 X 24 


30m 








7.3 


97 X 97 


30m 



where Ri {~ 9200 AU) is the radius where the temperature 
starts to be constant. 

Given these assumptions, the emissivity exponent mini- 
mizing the mass dispersion of the interferometric measure- 
ments shown in Fig.^ is y6 = 0.85 + 0.15, yielding /cumm - 
0.025+0.006 cm^ g At larger radii with the single-dish mea- 
surements, we get jS = 1.1+0.6, which is poorly constrained be- 
cause these measurements are not spread out enough in wave- 
length. In the following, we use /3 = 0.85. The best power-law 
fits to the masses derived from the interferometric and single- 
dish measurements are M(r) = 1.9 x (r/4200 AU)"''^ Mq and 
M(r) = 3.1x(r/4200AU)°'^^Mo,respectively(seeFig.n. Afit 
to an the data gives M(r) = 2.7x(r/4200AU)° " Mq. For com- 
parison, the mass of the singular isothermal sphere enclosed in 
4200 AU at 10 K is only 0.34 Mq. We note that the interfer- 
ometric measurements should resolve out any extended emis- 
sion from the environment, and therefore give a more reliable 
estimate of the mass distribution within 13 " in the IRAS 4A 
envelope. 

Assuming a spherical geometry and a power-law density 
profile - «o(^)^; the fit to the interferometric mass 
measurements yields p - 2.3 and riQ - 2.1 x 10^ cm"-^ at 
ro = 4200 A U. By analyzing their 2.7 mm i nterferometric 
measurements. iLoonev. Mundv. & Welch (2003) found no ev- 
idence for the presence of a circumstellar disk and they mea- 
sured a power-law index p between 1.7 and 2.3, which is con- 
sistent with our result (see their Tab. 1). Note also that no is 4 
times larger than the density of a singular isothermal sphere at 
10 K, and this ratio is even larger at smaller radii. 



by ISandell et alJ ill99ll) becomes 11.6 Lq. We use a dust 
emiss ivity a:^ = 0.1 x (A/250 fjmy^ cm^ g^' ( iHildebrandl 
Il983l) . The dust temperature profile is assumed to be fixed 
by the heating by the central protostar as suggested by 
iTerebev. Chandler & Andr'i ( ll993l) : 

with q - 2/(4 + /3) and L = 11.6 Lq. We compute the mass 
within a spherical radius R assuming optically thin emission 
and a temperature equal to Td{R). However, this calculation 
overestimates the true mass because of the radius dependence 
of the temperature. When the temperature and the density are 
single power laws with slopes -q and -p respectively and the 
Rayleigh- Jeans approximation is valid, the correction factor is 
(1 - 3^). For radii smaUer than ~ 9200 AU, Td(r) > 10 K and 



3.2. Spectroscopic signature of collapse 

In Fig. 13 we compare the CS, C^'^S, N2H+, HCO+, H'^CO-' 
and HC''^0^ spectra which we obtained toward the central po- 
sition of IRAS 4A with the IRAM 30m telescope, as well as 
CS, C^'^S, HCO+, H'^CO+ and DCO+ spectra which we took 
from the JCMT archive. Among the latter, CS(5-4) and CS(7- 
6) were observed at a position slightly off'set in right ascen- 
sion (less than 3"). We used the CS and C"^'*S frequencies mea- 
sured by iGottlieb et all (l2003h the N2H+(l-0) frequency de- 
termined b^ bore etalJ 12004 the DCO^ frequencies mea- 
sured by ' Caselli & Dorel ll2005l) . the HC O^ frequencies from 
the CDMS database (as of October 2005, Muller et alJf2005h . 
and the H'^CO'^ and HC '^O^ frequencies determined by 
ISchmid-Burgk et alJ (l2004 . 

The spectra of Fig. |2 show the classical signature of col- 
lapse, namely self-absorbed optically thick lines, with the 
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Fig. 2. Spectra toward the center of the IRAS 4A envelope obtained with the 30m telescope or taken from the JCMT archive: a) 
CS and N2H+ spectra, b) C^^s spectra along with CS(2-1) and N2H+(l-0), c) HCO+, H^^CO^, HC^^O+ and DCO+ spectra along 
with N2H^(l-0). The vertical dotted lines show the position of the blue peak, the red peak and the dip of the CS(2-1) spectrum 
in a) and b), and the same for HCO^(3-2) in c). The temperature scale is in Tmb- The beamwidth (HPBW) is given for each 
spectrum. 
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Fig. 3. a) CS(2-1) intensity maps towai'd IRAS 4A integrated 
over [-4.10,5.91] km s"' (thick contours), [6.06,8.31] km s"' 
(grey scale), and [8.46,15.50] km s"' (dashed contours). The 
contours go from 1.5 to 6 by 1.5 and from 1.5 to 3 by 1.5 
K km s"', respectively, b) CS(5-4) intensity maps integrated 
over [-1.37,5.91] km s"' (thick contours), [6.26,8.11] km s"' 
(grey scale), and [8.46,12.00] km s"' (dashed contours). The 
contours are 2 and 4, and 0.7 and 1.4 K km s respectively. 
The intensity scale is T*. The IRAM 30m beam size is shown 
in each panel (HPBW). The big/white and small/black crosses 
mark the positions of IRAS 4A and IRAS 4B respectively. 



blue peak being stronger than the red one, and low op- 
tical depth spectra pe aking in between (see lEvansI Il999t 
[Mvers, Evans & Ohashil EoOOl) . This signature is seen at 20 " 
in the East/West direction in our HCO^(3-2) map, but is only 
seen on the central position in our CS(5-4) map with a 10 " 
spacing. Since the smoothing of our CS(5-4) map to the JCMT 
resolution yields a spectrum very similar to the JCMT one, the 
infall signature seen in the JCMT spectrum comes entirely from 
the inner 10 " probed by the IRAM 30m telescope. 



In the c lassical spectroscopic signature of collapse motions 
llEvanslfl999<) . the shift of the dip of the optically thick lines 
with respect to the systemic velocity of the source given by the 
centroid velocity of an optically thin line is an indication of the 
infall velocity of the absorbing material. A fit to the hyperfine 
multiplet of N2H+(l-0) yields a low optical depth of 0.33 ±0.02 
for the isolated component 101-012, and a centroid velocity 
^N2H+ = 7.24 + 0.01 km.s"', which we take as the systemic 
velocity of the source. With respect to this systemic velocity, 
we measure shifts of 0.4, 0.4 and 0.5 km s"' for CS(2-1), (3- 
2) and (5-4), respectively, and shifts of 0.5 and 0.5 km s"' for 
HCO^(3-2) and (4-3), res pectively, all being measur ed with an 
uncertainty of 0. 1 km s~' . IPi Francesco et alJ ( I2OOII) measured 
a self-absorption at 7.64 km s~' with the Plateau de Bure in- 
terferometer in H2CO(3i2-2n), i.e. at the same velocity as the 
self-absorption in our CS(2-1) spectrum. 

The C"''*S(2-1) and (3-2) spectra are also asymmetric, with 
even a marked dip for C^^'*S(2-1). These lines could be optically 
thick and self-absorbed, like the CS lines. On the other hand, 
they could be optically thin but probe two velocity components 
along the line of sight. This would be the case for an infalling 
envelope with strong deplet ion at the center (see the case of 
LI 544 in I Caselli et alJl2002h . or with a small non-thermal ve- 
locity dispersion compared to the infall velocities, or for two 
independent velocity components, one being physically unre- 
lated to the envelope. We will try to address this question in 
Sect. |3and|5l With a signal-to-noise ratio of 3, only one peak 
is detected in the JCMT C"'"*S(5-4) spectrum, at the same ve- 
locity as the C-''*S(2-1) blue component. 

3.3. Extended wings tracing the outflow 

In Fig. 13 we present contour maps of the CS(2-1) and CS(5-4) 
line wing emission overlaid on the intensity maps of the line 
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Fig. 4. Map of H'^CO+(1-0) spectra taken toward IRAS 4A 
with the IRAM 30m telescope. Each spectrum is divided into 
two components, filled in black for v < 7.30 km s ' and in grey 
for V > 7.30 km s ' . The half-power beam width is 28". 

core. Both show a clear bipolar morphology whose direction 
RA. ~ 10° matches the direction of the outflow known in this 
source at small scale (see Sect. 0. The wings of the CS (and 
HCO^) central spectra in Fig.|2lare thus certainly associated 
with the molecular outflow. They extend to larger velocities in 
CS (from V = -2 to 15 km s ') than in HCO+ (3 to 1 1 km s '). 

3.4. Maps of low optical depth tracers 

Figure 13 shows a map of H'^CO^(l-O) spectra taken toward 
IRAS 4A with the IRAM 30m telescope. Nearly all spectra 
are double-peaked, as are also our C^^S spectra with a slightly 
worse signal-to-noise ratio. The intensity of the blueshifted 
component peaks at the IRAS 4A central position, this be- 
haviour being more pronounced in the North-East/South-West 
direction than in the perpendicular direction. This component 
is thus certainly physically associated with the IRAS 4A en- 
velope. By fitting two Gaussian components to each spectrum, 
we measure for the blueshifted component a centroid veloc- 
ity gradient of ~ 9.7 km s"' pc"' with a position angle RA. 
= 38° from North to East. On the other hand, the intensity of 
the redshifted component increases from the South-West to the 
North-East, and shows only a small velocity gradient of 2.4 
km s"' pc"' at P. A. = 7°. The interpretation of these velocity 
gradients in term of rotation is not straightforward. 

In Fig. |5lwe present large-scale maps of N2H^(101-012) 
intensity integrated over the same velocity ranges as for 
H'-'CO^(l-O) above (see Fig. 13 . First of all, we notice that 
the N2H^(I0I-0I2) emission has a local maximum very close 
to IRAS 4A in both Fig. [5^ and b, which strongly suggests 
that both the blueshifted and redshifted components are related 
to the protostar Second, these maps show that the N2H^(101- 
012) emission traces the sam e structu re as the dust continu um 
emission (e.g. ISandell & Knee..2001t iHatchell et allbOOSi) : a 
filamentary structure extending along a direction with a posi- 
tion angle RA. ~ -45°, which bends at RA. ~ -90° in the West. 
Furthermore these maps also provide additional kinematical in- 
formation and show that the two components which we iden- 
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Fig. 5. Maps of N2H+ (101-012) intensity integrated over a) the 
blueshifted velocity range [5.95,7.30] km s ' and b) the red- 
shifted velocity range [7.30,8.65] km s"', observed with the 
IRAM 30m telescope (HPBW 26 ", shown on the left). The 
first contour and contour step are 0.25 K km s ' in T* scale in 
both maps. The big cross marks the central position which we 
used for IRAS 4A. The small crosses mark the position o f the 
VLA a nd BIMA continuu m sou rces observed by Reipurth e t alJ 
ll2002i) and lLoonev et a iJ 12003) . The dashed line separates the 
two components. It emphasizes the spatial shift between them 
and their shell-like structure. 

tified in H'^CO^(l-O) have the same morphology but are spa- 
tially shifted: the "blue" filament is shifted by ~ 30" toward the 
South- West with respect to the "red" filament. Since they share 
the same morphology, they are likely to be physically related 
and not overlap by chance along the line of sight. 

4. Comparison to collapse models 

The classical signature of infall seen in our CS and HCO^ 
observations (see Sect. 13.21 and Fig. |3 suggests that the en- 
velope of IRAS 4A is collapsing, confirming the analysis of 
Di Franc esco et alJ (^001). The latter suggested that the in- 
verse P-Cygni profiles which they observed with the IRAM 
Plateau de Bure Interferometer toward IRAS 4A in H2CO(3i2- 
2ii) were created by an infalling envelope with a mass infall 
rate of 1.1 x lO"'* M© yr"'. This rate is nearly two orders of 
magnitude larger than the standard accretion rate at 10 K 
(.Shu. 1977 ). and is the result of both larger densities and veloc- 
ities (see Sect.lTland B.U . It is thus very unlikely that the col- 
lapse of this protostellar envelope has occured spontaneously, 
and we investigate in this section models of collapse induced 
by a fast external compression. 

4.1. Model of collapse induced by compression 

iHennebelle eTal] ll2003[ ! 2004) investigated with a smoothed- 
particle-hydrodynamics code the collapse of prestellar and pro- 
tostellar cores driven from the outside by an increase of the 
external pressure. In the following, we compare our observa- 
tions to their non-rotating model with a rapid compression, i.e. 
= 0.03, where <p is the ratio of the time-scale on which the ex- 
ternal pressure doubles to the initial sound-crossing time. Such 
a rapid compression is required to get a density enhancement 
matching the density profile of the IRAS 4A envelope (see 
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Fig. 6. a) Density, b) infall velocity and c) mass infall rate pro- 
files of the collapse model induced by compression described 
in Sect. Oat times 0.8402, 0.8533, 0.8581, 0.8657, 0.8763, 
0.8920 and 0.9249 xlO^ yr, corresponding to renormalized 
times since the formation of the central protostar of approxi- 
mately -7, 7, 11, 19, 30, 45 and 78 xlO^ yr, respectively. The 
thick line corresponds to 1.9 xlO'* yr and indicates the "best- 
fit" model used in Sect. 14.31 and shown in Fig.0 The dashed 
areas in (a) show the constraints deduced for IRAS 4A from 
the interferometric and single dish continuum measurements, 
respectively (see Sect. 13. 1> . 



Sect. 14. 2> . On the other hand, a more rapid compression pro- 
duces infall velocities too large to get a good agreement with 
our observations toward IRAS 4A. The initial conditions of the 
model correspond to a core in stable hydrostatic equilibrium. 
It is embedded in a hot rarefied external medium and has the 
structure of a Bonnor-Ebert sphere with a dimensionless radius 
= 3, i.e. smaller than the dimensio nless radius (i, - 6.45 of 
the critical Bonnor-Ebert sphere (e.g. iBonnoJI 1 956 l). We nor- 
malized the model to a kinetic temperature of 10 K and a total 
mass of 6.0 Mq. 

4.2. Mass distribution 

Figure |6t compares the evolution of the density profile of the 
model described in Sect. l4.1l with the density profiles deduced 
from the interferometric and single-dish continuum measure- 
ments toward IRAS 4A in Sect. 13. II The fast external compres- 
sion is responsible for a significant increase of the density pro- 
file with respect to the profile expected for a singular isother- 
mal sphere. The agreement with the profiles deduced from the 
observations is good in the range 1500-10000 AU at times 1-3 
xIO'* yr since point mass formation, which are of t he same or- 
der a s the estimated lifetime of Class protostars ( i Andre et alJ 
I2OOO1) . The disagreement at smaller radii, the model being less 
dense than what we deduced from the observations assuming 
spherical symmetry, could be due to the presence of an unre- 
solved disk not removed from the interferometric fluxes ana- 
lyzed in Sect. 13. II However, Loonev, Mundv, & Welch (2003) 
did not find evidence for the presence of such a disk. At larger 
radii, the single-dish measurements might be contaminated by 
large-scale emission from the cloud along the line of sight, 
not related to the protostellar envelope itself, although the dual 
beam technique used for the observations should have removed 
most of it. 



4.3. Radiative transfer modeling 

We used the Monte-Carlo-based radiative transfer code 
MAPYS O in one spherical dimension (cf. Blinder 1997; 
iBelloche et alJ2002h to compute the CS, C^'^S and N2H+ spec- 
tra that would be observed for a source set in collapse by a fast 
compression, as described by the model of Sect. 14.11 In addi- 
tion to the density and velocity profiles of the hydrodynamical 
model at time 1.9 xlO"* yr (see Fig.|6]and^, e), we used a 
uniform CS/C-^'^S isotopic ratio of 22.5 but non-uniform abun- 
dance profiles (see Fig. Et). For N^H ^, we used the collision 
rates of HCO^ with H2 llFIow e and assumed a constant 

abundance of 4.1 x 10 '°. We did not model the hyperfine struc- 
ture. We simply rescaled the abundance effectively used for the 
calculations by the statistical weight ratio ^ and fitted only 
the 101-012 optically thin component. For all molecules, we 
added a non-uniform non-thermal line broadening to the ther- 
mal broadening (see Fig. 0J)- Although the hydrodynamical 
simulation was isothermal, we used a non-uniform kinetic tem- 
perature to compute the spectra which are very sensitive to the 
heating by the central protostar The temperature was set to the 
dust temperature profile described in Sect. l3.1l (see Fig. 05). 

In Fig.|8lwe compare the synthetic spectra at time 1.9 xlO'* 
yr to the spectra observed toward IRAS 4A along the longi- 
tude axis, i.e. roughly perpendicular to the outflow axis (see 
Fig.EJ. The asymmetry and the absorption dips of the central 
CS(2-1) and CS(3-2) spectra are well reproduced. The asym- 
metry of these optically thick lines is enhanced toward the West 
and reduced toward the East at + 20 ", which is not accounted 
for by our ID spherical model. However including some ro- 
tation could improve the modeled spectra in that respect (see 
Belloche et al. 2002 for such an eff^ect in IRAM 04191). The 
line wings are not reproduced by our model, which was ex- 
pected since they trace the ouflow that we did not include in 
our modeling (see Sect. 13. 3> . The relative contribution of the 
outflow gets stronger with the upper energy level of the tran- 
sition, and even masks the red peak of the CS(5-4) spectrum. 
This spectrum is reasonably well matched by the model, al- 
though the shape of the dip is not fully reproduced. 

The central C^''S(2-1) modeled spectrum matches well the 
observed one. It is optically thin (r = 0.4) and the two peaks re- 
sult from the large infall velocities dominating the non-thermal 
line broadening in the inner parts of the envelope. We empha- 
size here the crucial importance of modeling optically thin lines 
of a less abundant isotopologue simultaneously with the op- 
tically thick lines of the main molecule. This indeed reduces 
the space of free parameters drastically - especially the non- 
thermal line broadening and the abundance - and helps to 
derive reliable infall velocities from the optically-thick, self- 
absorbed, asymmetric spectra. In particular, a first strong con- 
straint is set on the abundance profile by the spatial varia- 
tions of the C^'*S(2-1) integrated intensity, which imply here 
a decrease of the abundance toward the center. Besides, a sec- 
ond constraint on the abundance in the external parts comes 
from the depth of the absorption dips of the CS(2-1) and 
CS(3-2) lines. As a result, this set of data reveals that CS is 
strongly affected by depletion in the IRAS 4A enveloppe (see 
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Fig. 7. a) Density, b) kinetic temperature, c) CS molecular abundance, d) non-thermal line broadening (FWHM), and e) infall 
velocity profiles of the collapse model at time 1.9 x lO"* yr used for the radiative transfer modeling of Fig.|8l 
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Fig. 8. CS, C S and N2H+ spectra (histograms) in units of 
main beam temperature observed toward IRAS 4A along the 
longitude axis (see references in Fig.|2}. Synthetic spectra cor- 
responding to the collapse model at time 1.9 x 10"* yr described 
in Sect. 14. H are superimposed in red/thick line. The dotted line 
indicates the envelope systemic velocity used for the computa- 
tion of the radiative transfer (7.20 km s ' for CS and C^'^S, and 
7.24 km S-' for N2H+). 



iBelloch e et alJ2002l for a similar result in the Class protostar 
IRAM 04191). 

The N2H^(101-012) line is also well reproduced, although 
the observed blue peak is not as blueshifted as in the model. 
The opacity of the modeled spectrum is 0.24. It is in reason- 
able agreement with the opacity obtained from the Gaussian fit 
to the hyperfine structure (see Sect. 13. 2t . and the small discrep- 
ancy may be due to our simplified treatment of the radiative 
transfer of N2H^. With the uniform abundance, the N2H^(l-0) 
transition is more sensitive to the inner region than the C^'*S(2- 
1) hne. We tried a model including depletion of N2H^ above 
~ lO"' cm"^ (see lBelloche & Andre 2004 forNTH+ depletion in 



IRAM 04191) but the shape of the line was worse, the blue and 
red peaks being much more pronounced than they are in the ob- 
served spectrum. The reason why the uniform abundance pro- 
duces a profile with less pronounced blue and red peaks is that 
in this case, the beam peaks up more material at small velocities 
projected along the line of sight, i.e. the material infalling along 
a direction close to the plane of the sky. This shows that the 
double-peaked profile of the C^'*S(2-1) line is not only due to 
the large infall velocities compared to the non-thermal broad- 
ening, as stated above, but also to the depletion in the central re- 
gion. In our model, the shape of the N2H^(101-012) line there- 
fore indicates no (or very little) depletion of N2H^ in the inner 
regions. 

As we noted above, the non-thermal line broadening in the 
inner parts of the envelope is strongly constrained to be small 
(AVnt ~ 0.4 km s"', i.e. ctnt ~ 0.17 km.s"') by the double- 
peaked line shape of the low-optical-depth C^'*S(2-1) spectrum. 
The model with AVnt =1.1 km s"' in Fig.|9lshows indeed that 
a larger broadening smoothes out the C-''*S(2-1) double-peaked 
line profile which traces the infall velocities of the front and 
rear hemispheres. It also broadens the N2H^(101-012) spec- 
trum too much and produces wings which are not observed. 
On the other hand, the CS(2-1) and (3-2) broad absorption dips 
require a much larger non-thermal line broadening in the exter- 
nal parts of the envelope where they are produced (AVnt ~ 1 -2 
km s"', i.e. ctnx ~ 0.5 km.s"'). The model with AVnt - 0.4 
km s ' in Fig. |9l shows indeed that a smaller broadening can 
not produce the broad absorption dips observed in CS(2-1) and 
CS(3-2). We have therefore a strong evidence that the non- 
thermal line broadening is not uniform in the IRAS 4A en- 
velope, decreasing by nearly a factor of three toward the cen- 
ter, from supersonic down to subsonic values. This conclusion 
agrees quite well with the small linewidths AVnt ~ 0.5 k m s~^ 
observed toward IRAS 4A in N2H-i-(l-0) by Di Franc esco et alJ 
(2001) at small scales with PdBI (r ~ 1000 AU) and the 
much larger linewidths AVnt ~ 1.0-1.9 km s"' measured in 
the NGC 13 33 molecula r cloud with the low-density tracer 
C"*O(l-0) bv lWarin et al1lll99(il) . 

The time elapsed since point mass formation in the col- 
lapse model is well constrained by the peak separations and the 
Unewidths of the CS and C^'^S spectra. The infall velocities of 
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Fig. 9. CS, C^"^S and NaH^ spectra in units of main beam tem- 
perature observed toward the center of IRAS 4A (histograms). 
Each column shows a different model {thick line). Synthetic 
spectra corresponding to the collapse model at time 1 .9 x 10"* yr 
described in Sect. 14. H are superimposed in the central column 
{in red). In the left and right columns, the synthetic spectra {in 
green and blue) were obtained with the same input model ex- 
cept for the non-thermal broadening which was in both cases 
uniform, with AVnt = 0.4 and 1.1 km s"', respectively. The 
dotted line indicates the envelope systemic velocity used for 
the computation of the radiative transfer (see Fig.|8}. 



the model at time 7 x 10^ yr in Fig. [TO|are not large enough 
to produce the right separation of the peaks of the CS(3-2) and 
C"'^S(2-1) lines. On the other hand, the much larger infall ve- 
locities of the model at time 4.5 x 10"^ yr produce too broad 
spectra. As a result, we find a good match for a time of 1-2 
xlO^ yr elapsed since the formation of the central protostar. 

Finally, we could not find in the framework of our ID 
spherical modeling a better match for the C''^S(3-2) and 
C^'*S(5-4) central spectra. The modeled spectra of Fig. |S] are 
not asymmetric enough for these two lines. Since our model- 
ing suggests low optical depths of 0.35 and 0.04 for these two 
lines, it is very unlikely that the observed asymmetry results 
from an optical depth effect. This rather suggests that the ma- 
terial emitting in these two lines at the blue and redshifted ve- 
locities - i.e. in the rear and front hemispheres, respectively, in 
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Fig. 10. CS and C^'*S spectra in units of main beam tempera- 
ture observed toward the center of IRAS 4A {histograms). Each 
column shows a different model {thick line). Synthetic spec- 
tra corresponding to the collapse model described in Sect. 14.11 
at times -7, 7, 19 and 45 xlO^* yr are superimposed from left 
to right {in light blue, blue, red and green). For each model, 
the scaling of the abundance profile was adjusted to match the 
peak intensity of the C^'*S(2-1) line. The dotted line indicates 
the envelope systemic velocity used for the computation of the 
radiative transfer (see Fig.|S}. 



the framework of this collapse model - is denser for the former 
and less dense for the latter than assumed in the model. 



5. Discussion 

5.1. A vigorous collapse induced by a fast 
compression 

In Sect. 0] we found a very good agreement between the con- 
tinuum and CS/C^'*S/N2H+ line observations toward IRAS 4A 
and the collapse model of Hennebelle et al. (2003) induced by 
a fast external compression with (p - 0.03. The best match is 
obtained for a time of 1-2 xlO^ yr elapsed since the forma- 
tion of the central protostar. This corresponds to a mass in- 
fall rate of 4-7 xlO"^ Mq yr"' in the range of radii [10\lO'*] 
AU (see Fig.|5J;). However the actual mass infall rate is proba- 
bly even larger since the density profile of our model does not 
match perfectly (and underestimates by about 25%) the den- 
sity profile deduced from the observations (see Fig.|6t). Note 
also that the dust opacity we used in Sect. 13.11 is 2.5 times 
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larger than the opacity usually a ssumed for protostellar en- 
velopes ("e.g. lMotte & Andrel200ll) . Therefore we estimate that 
the actual mass infall rate is most probably in the range 0.7-2 
xlO"^ Mq yr This is in very good agree ment with the mass 
infall r ate of 1 . 1 x 10 yr"' deduced bv lPi Francesco et alJ 
il200lh from their PdBI observations. As we already mentioned 
in Sect. 0] such a mass infall rate of 1 xlO""^ Mq yr"' is 60 
times larger than the standard accretion rate ^ at 10 K jShj 
Il977l) . It is the result of both higher densities and velocities pro- 
duced by the fast compression wav e propagating inwards used 
in our model. IFatuzzo et alJ (l2004l) found indeed that the mass 
infall rate scales roughly linearily with the initial overdensity 
and the initial velocity field in their self-similar collapse solu- 
tions where the dense core is initially either overdense com- 
pared to the singular isothermal sphere or has nonzero initial 
velocities. 

In the framework of the self-similar inside-out collapse 
model with an effective sound speed including a non-thermal 

contribution, the "standard" accretion rate would be ' ^'^ — 
(Shu et al. 1987). Using the non-thermal velocity dispersion 
(Tnt = 0.17 km s ' derived in the inner parts of the enve- 
lope in Sect. 14.31 which is typical of the level of turbulence 
in prestellar condensations (e.g. Belloc he et al. 2001) but prob- 
ably even an upper h mit for IRAS 4A since the smallest dis- 
persion measured by ' Pi Francesco et all (12001 ^ with PdBI is 
0.09 km s the "effective" accretion rate would be at most 3 
times larger than the pure isothermal one and cannot account 
for the large mass infall rate measured in the IRAS 4A enve- 
lope. Ori the other hand, the L arson-Penston similarity solution 
jLarsonl 1 969t IPenstonI 1 969l) has a uniform mass infall rate of 
29 X at point mass formation and its extension after point 
mass formation by Hunter ( 1977i) has a mass infall rate ap- 

proaching asymptotically 47 x This value is closer to our 
result for IRAS 4A, which is not surprising since among the 
numerous similarity solutions of isothermal gravitational col- 
lapse the Larson-Penston solution can be physically interpreted 
as a model associated with a st rong external compression wave 
dWhitworth & Summerslfl98^ . However, although helpful to 
understand the physics of gravitational collapse, similarity so- 
lutions are not well suited for detailed comparison with obser- 
vations since they do not have realistic boundary conditions. 

By contrast to the large mass infall rate measured in the 
IRAS 4A envelope, the spontaneous collapse undergone by 
the very young Class protostar IRAM 04191 in the Taurus 
molecular cloud is much less vigorous with a mass infall rate 
of 2 X ^ only (see iBelloche et aljEoo3) . This could explain 
the order of magnitude difference between the bolometric lu- 
minosities of the two protostars which have yet approximately 
the same age since the b eginning of the main accretion phase 
rseelT,esaffreetal.l2005l). 

5.2. Origin of the external compression 

The N2H^(l-0) maps of Fig.|5]show that the material emitting 
at "blueshifted" velocities is displaced toward the South- West 
with respect to the material emitting at "redshifted" velocities 



and that both components have the same shel l -like m orphol- 
ogy, at the edge of cavity 2 of iLefloc h et alj (| l998h seen in 
maps of continuum emission (see also lSandell & KnedboOlt 
iHatchell et al. 2005). We interprete the velocity difference be- 
tween the two components as the result of "inward" motions 
and conclude that the "blueshifted" material is in the back- 
ground and moving toward the "redshifted" material along a 
direction projected onto the plane of the sky at PA. ~ 45° in 
the IRAS 4 region and PA. ~ 0° in the SVS13/VLA 12 re- 
gion. These geometry and velocity structure remind us of an 
expanding shell, the center of which would be located in the 
background and South- West of NGC 1333. The "blueshifted" 
component would be the expanding shell itself while the "red- 
shifted" component would be the ambient medium. In this sce- 
nario, the whole region around IRAS 4 and IRAS 2 should 
show signatures of "inward" motion, i.e. the motion of the 
"blueshifted" component toward the "redshifted" component. 
This might be indeed the case since 'Wals h. Bourke & MversI 
(2006) claim that they have detected "global infall" in this re- 
gion. On smaller scales in the IRAS 4A region, this scenario fits 
also, from a kinematical point of view, into the framework of 
the collapse model we investigated in Sect.|4] The external per- 
turbation triggering the collapse of the dense core would then 
be the impact of this expanding shell on the ambient medium. 

In this scenario, the origin of the external perturbation 
should be located in the background and South- West of the 
IRAS 4/SVS13 region. The perturbation could result from the 
expansion of an H II region, from a stellar wind or from a pro- 
tostellar outflow. A nearby bright star in the South- West is re- 
quired for the first two assumptions. The star BD +30°547, with 
a spect ral type between A 6 V and B 7 V ( Aspin 2003, although 
ICernisI 1990 f ound a type G 2 IV), would be the best candi- 
date. However. lPreibischl ll2003h found with XMM-Newton that 
this star is in the foreground, so it cannot be, in our scenario, 
the origin of the external perturbation in the IRAS 4 region. 
Since NGC 1333 is filled with many outflows ("k nee & Sa ndelll 
2000), the coll apse of IRAS 4 could have been triggered b y one 
of them. Like lHennebeUe et alJ ll2003h . iMotovama & Yoshidal 
(2003) investigated numerically the collapse of dense cores 
triggered by an external perturbation. They initiated the col- 
lapse by setting a nonzero inward velocity at the boundary of 
the dense core. They found that the maximum value of the mass 
infall rate is proportional to the momentum given to the dense 
core by the external perturbation, and in particular that an in- 
put momentum of ~ 0.1 Mq km s"' yields a mass infall rate 
of ~ 1 X lO"'* Mq yr"'. IRAS 4 A could thus have been set in 
collapse by an external perturbation with a momentum of ~ 0.1 
M0 km s"' . Such a momentum is a typical value for the present 
molecular outflows in NGC 1333 (Knee & Sandefl 2000). The 
collapse of IRAS 4A could therefore have been triggered, di- 
rectly or indirectly, by the shock created by a present or former 
outflow. This would fit into the scenar io of shock driven se- 
quential star formation put forward by IWarin et'aP ( 1 19961) in 
NGC 1333. Indeed, the cavity in the Sou th- West of IRAS 4A 
could be the relic of a former outflow (see lOuillen et alJ2005l) . 
However, we have not yet found a good protostellar candidate 
which could have been the origin of such an outflow in the 
South-West. 
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We mentioned the limitation of our ID spherical model- 
ing for IRAS 4A in Sect. 14.31 and suggested an asymmetry of 
the source along the line of sight, the rear hemisphere being 
denser than the front hemisphere. This asymmetry could result 
from the fast external compression itself, which i s unlikely to 
be isotropic at the scale of the IRAS 4A envelope. IBossI i ll 9951) 
showed examples of such asymmetries in dense cores set in col- 
lapse by a shock wave hitting only one hemisphere (see his Fig. 
2). It would be geometrically consistent with the idea that the 
compression wave which set the IRAS 4A envelope in collapse 
is coming from the back, as deduced above from the morphol- 
ogy and kinematics of the N2H^(l-0) maps. 



5.3. Non uniform turbulence 

We found in Sect. 14.31 that the non-thermal line broadening, 
which is likely to result from turbulent motions, decreases to- 
ward the center of IRAS 4A from supersonic to subsonic ve- 
locities by nearly a factor of 3. A similar conclusion was al- 
ready drawn by Di Francesco et al. ( 2001 ) based on their N2H^ 
interferometric data. The turbulence is thus supersonic in the 
outer parts of the IRAS 4A envelope but only subsonic in the 
inner dense region. The subsonic turbulence in the inner part 
of the envelope is reminiscent of the low level of turbulence 
found in prestellar condensations (e.g. in the p Oph protoclus- 
ter, see [Belloche et al. 2001). It could then indicate the con- 
ditions prevailing in the dense core before the onset of col- 
lapse. In the outer parts of the envelope, the supersonic tur- 
bulence (Tnx ~ 0.50 km s"' which we derived from the CS 
self-absorption is nearly as l arge as the a mplitude of the non- 
thermal motions deduced by IWarin etaP ( Il996) from the low 
de nsity tracer C^^O(l- O) in the molecular cloud (see also Fig. 7 
of lOuillen et'ai]l2005l for C'^O data with a better spatial res- 
olution). The outer parts of the IRAS 4A envelope are there- 
fore affected by the supersonic turbulence which permeates the 
whole cloud and which was ma ybe powered by th e numerous 
outflows located in NGC 1333 ( Ouillen etalJ2005l) . In the sce- 
nario we proposed in Sect. 15. 21 the supersonic turbulence could 
have been generated in the outer parts of the IRAS 4A envelope 
by the expansion of the south-western cavity. 



5.4. Rotation and binary formation 

Detecting rotation and measuring its magnitude in the envelope 
of IRAS 4A is not straightforward. We mentioned in Sect. 13.41 
that the blue component of the low optical depth tracers shows 
a velocity gradient of ~ 9.7 km s"' pc"' at RA. ~ 38° over 
~ 40", which is of the same order as the gra dient measured 
in IRA M 04191 that we interpreted as rotation ( Bel loche et alJ 
However, there are two caveats for IRAS 4A. First, the 
double-peaked structure of our single-dish low-optical-depth 
spectra makes the a n alysis in terms of rotation difficult. Second, 
iDi Francesco et alJ (l200ll) measured with PdBI in N2H+(1- 
0) centroid velocity differences of about 0.5-1 km s ' over 
~ 12", which corresponds to a velocity gradient of ~ 30-50 
km s ' pc ', at a position angle RA. ~ 135°. Their position 
angle was only a rough estimate and a fit to their centroid ve- 



locity map would actually give a position angle probably close 
to 90° (see their Fig. 2b), but this is still ~ 50° away from our 
single-dish measurement and it is unclear whether these veloc- 
ity gradients trace indeed rotation. 

However, if we assume that the velocity gradient we mea- 
sured in Sect. l3.4l is a good estimate of the amount of rotation 
in the envelope and suppose at first order that it is a solid-body 
rotation, then we find a ratio of rotation energy over gravita- 
tional energy of /? ~ 0.02. T his is similar to the conditions 
investigated bv iHennebelle et al . (2004) to study the fragmen- 
tation process in a rotating core set in collapse by a fast exter- 
nal compression. Their main conclusion was that a fast com- 
pression promotes fragmentation and the formation of multi- 
ple protostars. With a very rapid compression - 0.1), they 
showed that this process occurs through the formation of a 
ring which becomes very quickly unstable. It breaks up into 
several pieces which, in the example they showed, merge into 
two protostars. For their simulation with a dense core of 1 M©, 
the proto-binary separation is about 140 AU (see their Fig. 7), 
which corresponds to 840 AU if we normalize to the mass of 
the IRAS 4A envelope (6 Mq). IRAS 4A is indeed a binary sys- 
tem wit h a separation of 1.8", i.e. 570 AU at a distance of 31 8 
pc (see lLoonev. Mundv. & WeTchlbOOOt iReipurfli et al.ll2002h . 
It is therefore tempting to conclude that the formation of the 
IRAS 4A binary system results from the fast external compres- 
sion which has set up the dense core into collapse. 

5.5. Collapse of IRAS 4B 

Like toward IRAS 4A. lDi Francesco etaP ( l200ll) found an in- 
verse P Cygni profile toward the Class protostar IRAS 4B 
(Sandell et al. 1991). They concluded that it is also collapsing, 
with a similar large mass infall rate of 3.7 x 10"^ M© yr"'. 
Since it lies along the same N2H+ filament as IRAS 4A (see 
Fig. |5} and is in a similar evolutionary stage, we suggest that 
the gravitational collapse of its envelope was triggered by the 
same external perturbation, the expanding shell suggested in 
Sect. 15:^ 

We note that IRAS 4B is an 11" bina ry, and possibly 
a system with even higher multiplicity (see iLav et alJll995t 
lEoonev et al. 2000). However the main comp anion was not 
detected in continuum emission at 3.6cm by iReipurth et alJ 
('2002') and is probably much more evolved than the primary 
(Looney et al. 2000). This raises the question whether it really 
belongs to the IRAS 4B system or is physically unrelated and 
appears just by chance projection close to the Class protostar. 

5.6. Alternative interpretations ? 

In the past, it has been claimed that the NGC 1333 molec- 
ular cloud includes several independant v elocity components 



seen in emission and/or absorption (see Lang er et alJ 
■Choi. .200k jChoi et al...2004) . In p articular, QioLeLal 



rejected the infall interpretation of iDi Francesco et~ 



199t, 



( 2004 1 

2m 



and proposed instead a foreground absorbing layer at 8 km s ' 
physically unrelated to IRAS 4A. They assumed a very low 
systemic velocity of 6.7 km s"' for IRAS 4A, yielding highly 
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redshifted absorption dips in the optically thick spectra, and ar- 
gued t hat standard collapse models like the inside-out collapse 
model would never produce such highly redshifted 

dips given the age of the protostar, and that the faster Larson- 
Penston model would produce much too broad linewidths in or- 
der to fit these highly redshifted dips. However, the model pre- 
sented in Sect. |4] shows that a systemic velocity of 7.2 km s"' 
and a collapse induced by a fast external compression can fit 
both the optically thick and thin lines reasonably well. The 
fit to our data set does not require any additional absorbing 
layer p hysically unrelated to IRAS 4A as claimed bv lChoi et alJ 
(Hop). Our scenario, which requires only one source, is sim- 
pler. Therefore we think it is more likely. 

6. Summary and conclusions 

We have carried out a detailed analysis of the physical struc- 
ture of the envelope of the Class protostar IRAS 4A in the 
NGC 1333 molecular cloud. Our main results and conclusions 
are as follows: 

1. The density profile deduced from previous interferomet- 
ric and single-dish millimeter continuum measurements is 
steep (p = 2.3) and the envelope is much denser than a sin- 
gular isothermal sphere at 10 K. 

2. Our new set of molecular line data obtained with the 
IRAM 30m telescope shows the classical spectroscopic sig- 
nature of collapse with absorption dips of optically thick 
lines redshifted by 0.5 km s"' with respect to the systemic 
velocity of the source traced by optically thin lines. 

3. Both the density structure deduced from the continuum ob- 
servations and our set of CS, C^^S and N2H+ spectra can be 
well fitted by a model of collapse induced by a fast external 
compression with a time elapsed since point mass forma- 
tion of 1-2 X 10'* yr. This age is consistent with the lifetime 
estimated for Class protostars. 

4. We deduce a very large mass infall rate of 0.7-2 x 10"^ 
Mq yr ' in the envelope of IRAS 4A. It results from both 
the overdensity and the large velocities generated by the 
fast external compression. 

5. Our radiative transfer modeling reveals a strong decrease 
of the level of turbulence from supersonic velocities in the 
outer regions to subsonic velocities near the center of the 
IRAS 4A envelope. The subsonic turbulence is probably a 
relic of the conditions prevailing in the dense core before 
the onset of collapse. 

6. We found indications of the presence of an expanding shell 
in the vicinity of IRAS 4A which could be the origin of the 
compression. This expanding shell could have triggered the 
collapse of the nearby Class protostar IRAS 4B simulta- 
neously. However, we did not find evidence for the pres- 
ence of a nearby bright star which could have initiated this 
compression through the expansion of an H II region or a 
strong wind. On the other hand, the collapse of IRAS 4A 
and IRAS 4B could have been triggered, directly or indi- 
rectly, by an outflow, but we have not found a good proto- 
stellar candidate yet. 
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